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ABSTRACT 


Remote sensing technology has been widely applied in various fields, 
including oil, gas, and mineral exploration, spatial planning, and 
environmental monitoring. This paper describes the application of remote 
sensing technology for the potential study of a renewable micro hydropower 
plant (MHP) using Landsat 8 satellite data. The Sukaati Watershed, West 
Java, Indonesia, was selected as the case study area. Landsat 8 satellite data, 
acquired on August 21, 2020, was applied to extract information on land use, 
geology, and potential landslides. Drainage patterns, watershed boundaries, 
and head height were obtained from topographic map data. Drainage patterns, 
watershed boundaries, and land use are used to calculate flow rates. 
Geological map and landslide are the basis of layout of MHP components, 
such as water intake, dam, waterway, settling tank, penstock, and powerhouse. 
A field survey to acquire actual flow rate and head height was conducted to 
validate the results of the remote sensing data interpretation. Two potential 
sites of MHP were selected with a hydropower design of 129 kW and 5.18 
MW. This study showed that remote sensing technology is beneficial for 


studying the potential of MHP because fieldwork can be done more quickly 
and efficiently. 
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1. INTRODUCTION 

The level of electricity expenditure can be observed as an assessment measure of a country’s 
development level. The rise in the demand for energy worldwide is an affirmation that energy is the engine of 
growth in this 21st century [1], [2]. Renewable hydropower is an alternative source of energy that can meet 
energy needs and reduce dependence on fossil energy sources. Renewable hydropower is a clean, reliable, 
versatile, and low-cost source of electricity generation [2]-[10]. 

Hydropower potential studies carry huge uncertainties about the technical, economic, and 
environmental feasibility. Building a hydropower plant requires thorough studies of water drainage patterns, 
geological lithology, surface and subsurface structures. Image processing technology can play a vital role in 
the studies of the categories mentioned. Detailed field studies can be followed, once image processing 
accomplishes the initial stage of the studies. Remote sensing technology can generate, maintain, and renew 
up-to-date data of the earth's surface and its combination with socio-economic information. Remote sensing 
and GIS technology are very helpful for the hydropower plant feasibility study [11]-[16]. 
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Remote sensing and geographical information systems, with the development of satellite 
technologies, becomes possible to assess rapid and economic solutions for hydrology mapping, to determine 
a practical rainfall-runoff relation, water flow, particularly poorly gauged or ungauged basins [14], [17]-[28]. 
Remote sensing data is also useful for geological mapping. This method can minimize the cost and reduce the 
time of geological field surveys. The field survey can be done more effectively and more precisely [12], [13], 
[19], [21], [25], [29]-[33]. 

This study presents a methodology for utilizing Landsat 8 satellite data to study the potential for 
MHP. ER Mapper 7.3 and ArcGIS 10.3 are used for digital image processing and spatial data analysis. 
Utilizing satellite data through interpretation and analysis, geological, environmental (land use), 
hydrological, and other data are obtained to determine potential locations for MHP construction. The MHP 
study in the Cibuni river in Sukaati village is based on an estimate of the potential for electrical energy 
generated, which is quite large, which can be used for the tea processing industry, for the surrounding 
community who really need electricity for lighting and can be distributed to PLN (Indonesian state electricity 
company). The results of this interpretation are integrated with ground truth checks to validate to improve the 
accuracy of their interpretation further. Thus remote sensing data is proposed for the study of the potential of 
MHP in terms of minimizing costs and reducing the time of the preliminary survey of MHP. This paper also 
discusses the role of remote sensing in interpreting geology, land use, hydrology, and geological disasters, 
such as landslides in relation to the potential site of MHP, so that the sustainability of MHP infrastructure 
will be better maintained. 


2. RESEARCH METHOD 

Sukaati watershed, Cibuni river, West Java, Indonesia was selected as a case study area as shown in 
Figure 1. The Cibuni river originates on the western slopes of Mount Patuha II in the Rancabali district, 
Bandung regency, and empties into the Indian Ocean in Sindangbarang district, Cianjur regency. Ci Buni 
river has a length of 109 km and an area of 1,434.70 km2 covering three districts, the Cianjur, Bandung, and 
Sukabumi regency. With a large catchment area and a large discharge, the Cibuni River is a potential water 
resource that is used for various needs, including drinking water, irrigation, and flushing as well as 
generating electricity. 

Cibuni river flows along the southwestern region of the island of Java, which has a tropical 
rainforest climate (code: Af according to the Köppen-Geiger climate classification). The annual average 
temperature is around 18°C. The warmest month is October at an ordinary 20°C, and the coldest month is 
August at a regular 18°C. The average highest rainfall is in December at 570 mm, and the lowest is in 
September at 53 mm. The annual rainfall is high, between 2239 mm to 5779 mm. The average annual rainfall 
is 3766 mm. The preliminary reading of the RBI contour map showing drainage patterns, watershed 
boundaries, and potential heads identified two potential sites of MHP, Site 1 (Indragiri site), and Site 2 
(Sukaati site). 
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Figure 1. Study area, Site 1 (Indragiri site), and Site 2 (Sukaati site), Sukaati watershed, Cibuni River, West 
Java, Indonesia 
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By deploying remote sensing data interpretation such as mapping drainage pattern, watershed 
boundary, land use, geology, and potential landslide, several tasks of MHP potential studies can be covered. 
The research method is presented in Figure 2. The utilization of Landsat 8 satellite data was to create land 
use, geology, and potential landslide maps. The base map published by Geospatial Information Agency of 
Indonesia, with a scale of 1: 25.000, was the basis for making a topographic, drainage pattern, and watershed 
boundary map. 

The estimated flow rate was calculated from the drainage pattern, watershed boundary and land use 
spatial data. The head was estimated from the topographic map. Afterward, the geological and potential 
landslide map was used for layout design of MHP components, such as water intake, dam, waterway, settling 
tank, penstock, and powerhouse. Field measurement is performed to measure the head and water flow rate. 
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Figure 2. Research method 


The hydrological study aims to determine whether the available water discharge and plunge height 
can drive the turbine according to the desired power. The hydropower potential study includes measuring the 
minimum discharge flowing in a water channel/river, water discharge during a flood, and available waterfall 
height (height difference). The generated hydropower is calculated using: 


P=gxQxHxMn (1) 


where: P=estimated hydropower (kW), g=gravity (m/s*), Q=flow rate (m°/s), H=effective head (m), 
n=system efficiency=80 (%). 


2.2. Head and flow rate estimation 

Related to the head and flow rate calculation, the head, drainage pattern and watershed boundary 
can easily be identified from the digital elevation model (DEM) extracted from a topographic map. The 
vector data of the drainage patterns map can be overlaid to the raster image for comprehensive reading 
[34]-[38]. 
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Topographic studies include collecting data and information on the best location to obtain adequate 
water drop height, intake penstock, and tailrace. Based on topographic map, it can be identified two potential 
MHP site, Site 1, Indragiri watershed with a potential head of 25 m, and Site 2 Sukaati watershed with a 
potential head of 125 m. Parameters for calculating water flow consist of precipitation, catchment area, and 
type of land cover in the catchment area. Water flow is obtained from as shown in (2), (3): 


Q=0.00278 Cws x I x A (2) 
u Dici Ci Aj 
Cws = rT" Ai (3) 


where: Q is the water flow (m3/s), Cws is the runoff coefficient (depend on the land cover/land use), I is the 
precipitation (mm/hour), A is the watershed area (ha), Ai is the land area with land cover type 1, Ci is the 
runoff coefficient of the ground cover type 1, n is the number of land cover types. 

The precipitation data was provided by the rainfall gauge station data belonging to PTPN VIII Tea 
Plantation. Rainfall data is the average monthly rainfall obtained from measurement data for 20 years, from 
2000 to 2019. Land use is divided into four classes with a runoff coefficient values listed in Table 1. 


Table 1. Coefficient runoff in study area 


Land use/Land cover Runoff coefficient (C) 
Forest Land 0.10—0.40 
Agriculture Land 0.40-0.50 
Built-Up Land 0.50—0.70 
Bare Land 0.70—0.95 


The land use map was interpreted based on the land cover identification of remote sensing data. 
Land cover refers to the physical and biological covering up the land's surface, including water, vegetation, 
bare land, or artificial structures. Land use relates to human activities such as agriculture, forestry, as well as 
build up area. The land cover may be observed directly in the field or by remote sensing. Observations of 
landuse need to integrating social and natural scientific approaches to decide which human activities occur in 
various landscape parts, even when land cover appears to be the same. 

A Landsat 8 satellite data of the study area was downloaded from earth explorer-United States 
geological survey (USGS) with path type WRS 122 and Row 065. The acquisition date is August 21, 2020. 
The projection is UTM WGS84 Zone 46. Landsat 8 is the newest Landsat satellite launched and brings the 
operational land imager (OLI) and the thermal infrared sensor (TIRS) instrument. There are nine spectral 
bands of OLI and two spectral bands of TIRS. The Landsat 8 contains invisible multispectral bands and 
infrared bands. The layer of the multispectral bands consists of near-infrared (NIR), short-infrared (SWIR), 
red and green are employed as shown in Table 2. 


Table 2. Landsat 8 OLI/TIRS [39] 


Band Designations Wave Length Spatial Resolution 
Band 1 (Coastal Aerosol) (0.43-0.45 um) 30m 
Band 2 (Blue) (0.45-0.51 um) 30m 
Band 3 (Green) (0.53-0.59 um) 30m 
Band 4 (Red) (0.64-0.67 um) 30m 
Band 5 (Infrared) (0.85-0.88 um) 30m 
Band 6 (Short wave infrared) (1.57-1.65 um) 30m 
Band 7 (Short wave infrared) (2.11-2.29 um) 30m 
Band 8 (Pancromatic) (0.50-0.68 um) 15 m 
Band 9 (Cirrus) (1.36-1.38 um) 30m 
Band 10 (TIRS) (10.6-11.19 um) 100 m 
Band 11 (TIRS) (11.5-12.51 um) 100 m 


In this study, Landsat 8 OLI data is used for land use, geology, and landslides mapping. Land use 
mapping is carried out using a combination of band 5 (near infra-red/NIR), band 4 (Red), and band 3 (Visible 
Red). The Landsat 8 OLI was extracted to enhance the land cover features. The combination of bands can 
interpret the characteristics of the object, such as vegetation, water, and developed areas [40]. The 
normalized difference vegetation index (NDVI) represents the vegetation. Vegetation is a mathematical 
combination between the red and NIR band, which indicates vegetation's presence and condition. The 
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normalized difference water index (NDWI) represents the water by estimate the leaf water content and 
monitoring changes related to water content in water bodies. The NDWI uses a near-infrared and green band 
to enhance characteristics while eliminating soil and terrestrial vegetation. The normalized difference built- 
up index (NDBI) represents the built-up area. The NDBI uses a combination of the NIR band and SWIR 
band to highlight built-up areas automatically. The NDVI, NDWI, and NDBI were determined by utilizing 
the formulas (4), (5), and (6): 


NDVI=(NIR-Red)/(NIR+Red)=(Band 5-Band 4)/(Band 5+Band 4) (4) 
NDWI=(Green-NIR)/(Green+NIR)=(Band 2-Band 5)/(Band 2+Band 5) (5) 
NDBI=(Green-SWIR)/(Green+S WIR)=(Band 5-Band 6)/(Band 5+Band 6) (6) 


2.3. Geological study 

The geological study aims to convince various parties that the micro-hydro development program is 
technically feasible and suitable to support the development of MHP. Data and information required in the 
geological study include topographical and geological that will be used for the main building and drainage 
routes. The geological condition must support the quality of the core building, which consists of a dam, 
intake, water way, settling tank, penstock and power house. According to landslide prevalence, geology is a 
prime factor, hence the necessity to create a structural pattern map, faults map, and river channel map in a 
highly rugged volcanic mountainous terrain [41]. Geological data is obtained from a combination of image 
interpretation, topographic map interpretation and field checks to determine accuracy in the analysis of rock 
types and geological structures. Regional geological data refers to the geological map sheet Sindang Barang 
and Bandarwaru [42]. 

Remote sensing data is helpful for geological mapping [36], [43], [44]. In geological mapping, 
reflectance information of rocks in the SWIR and emissivity of the rocks in the TIR are noteworthy. In 
Landsat 8, SWIR is available in Band 6 and 7; TIR is available in Band 10, 11. In general, combinations 
involving bands from each spectral region (visible, mid-infrared, and SWIR II) display the most contrast on 
lithological that makes Landsat 8 (R=4, G=5, B=7) were enhanced further using histogram equalization 
contrast enhancement accordingly. In this study, geological and landslide mapping was carried out using a 
combination of band 4 (Red), band 5 (NIR), band 7 (Shortwave Infra-Red). Band rationing and combinations 
with most contrast for band ratios involving geology is a higher band divided by a lower one. Combinations 
(4/3, 5/2, 7/3) were selected for geological mapping. Field geological checks are carried out to determine 
rock types and other geological information related to micro hydro development. 


3. RESULTS AND DISCUSSION 
3.1. Flow rate estimation 

Landsat 8 OLI imagery in the form of a combination of visible wave and NIR bands is applied to 
obtain information on land cover and land use. Afterward, a vegetation index, water index, and build-up 
index are formulated. NDVI, NDWI, and NDBI data are displayed as a composite image and serve as a 
classification reference for making land use maps as shown in Figure 3 (a) and (b). An enhanced image 
provides better information for interpreting Landuse, Land cover as well as geology and landslide mapping. 

To obtain lithology information, Landsat 8 OLI imagery was used in the form of a combination of 
NIR and SWIR bands. Band rationing and assimilations (4/3, 5/2, 7/3) have the best contrast and form the 
input for the classification of geological mapping as shown in Figure 3 (c) and (d). 

The interpretation of Land-use is displayed in Land use map as shown in Figure 4. It appears that 
the study area is mostly plantation and forest, which has a low run-off coefficient. Land cover in the form of 
forests and vegetation causes a small fluctuation of flow rate, which means proper for the sustainability of 
MHP. Based on the remote sensing interpretation, the Indragiri watershed with 1,595 ha catchment area has a 
potential discharge of 1.30 m?/s, while Sukaati Watershed has a water catchment area of 6,515 ha and yields 
a 10.2 m°/s flow rate. The land cover in the study area is relatively homogeneous, namely in the form of 
forests in the upstream and tea plantations in the hills. The owner of the case study area is Forestry and Tea 
Plantation Company. Land tenure by Forestry and Tea Plantation Company will facilitate the management of 
water infiltration in the catchment area. Forest and tea plantations land cover will provide benefits to 
maintain the stability of the flow rate required for MHP, reduces the risk of erosion and landslides. 
Landslides mainly occur in steep river valleys due to vertical erosion of water. 
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3.2. Geological study and head estimation 

The physiography of the research area is the volcanic mountains of south Bandung, this is shown by 
its morphology such as cones and volcanic craters, then the shape of steep river grooves and valleys with 
heavy water flows. The geological interpretation is presented in geological map in Figure 5. The main river 
of the study area is the Cibuni River, which originates at Mount Patuha and Mount Kendeng, flows to the 
southwest, and empties into the Indonesian Ocean. Brook form steep and sheer incisions forming a radial 
pattern upstream and developing dendritic downstream. In the study area, all pyroclastic rocks are the 
structure of products of volcanic activity, namely Mount Kendang (QIk) in the northern part, which is an old 
volcano with cranial remains on it. Mount Patuha (Qvp) is on the eastern side with its parasitic crater in the 
form of a white crater and a lake of Patenggang is in the middle. The rocks consist of andesite lava, breccia, 
and sandy tuff containing pumice and lava deposits in Figure 6. River erosion tends to be vertical, indicated 
by the steep and sheer river valleys. The geological conditions of the study area provide an overview of the 
carrying capacity of the land for MHP development. It is of concern for civil construction design. The 
thickness of weathered soil and the high fertility rate of volcanic rocks make this area very suitable for tea 
plantations and community fields. 

The results of the Landslide potential interpretation is presented in landslide map in Figure 7. There 
are several locations with high landslide potential (the areas in red). Feature of potential Landslide in the 
image has a specific characteristic of the landslide crown, steep walls of the topographic valley, and rivers 
that differ from the surrounding slopes. Steep river valleys in soft and unstable rock types are a potential 
landslide zone. In volcanic rocks, generally, the weathering is very thick due to high rainfall and high slope 
levels. In general, soil movements occur in valleys or slopes composed of weathering soil from breccia and 
tuffs. In this study, areas with high landslide potential are generally not in the Sukaati catchment area. This 
result signifies that the Sukaati Watershed is the potential for the construction of MHP. Areas prone to 
landslides can be avoided by designing a more detailed site plant design. 





(c) (d) 


Figure 3. Show, (a) original image of Landsat 8 OLI (R=5 G=4, B=3), (b) enhanced image of Landsat 8 OLI 
(R=NDVI=(5-4)/5+4), GENDWI=(2-5)/(2+5), BENDBI=(5-6)/(5+6)), (c) original image of Landsat 8 OLI 
(R=4, G=5, B=7), (d) enhanced image of of Landsat 8 OLI (R=4/3, G=5/2, B=7/3) 
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Figure 6. Photo of outcrop on the cibuni river in the Sukaati area with steep river valleys, (a) andesite lava, 
(b) lava breccia 
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Figure 7. Landslide map of Sukaati watershed 


3.3. Field measurement results and comparative analysis 

Field measurement is performed to verify the results of remote sensing interpretation. Measurements 
were conducted in August 2020 (the dry season when discharge is minimum). The flow rate is obtained by 
multiplying the water velocity measured using a flow meter, and the cross-sectional area of the river. All the 
head and flow rate data from measurements and the remote sensing interpretation are then used to calculate 
the hydropower output as shown in Table 3. 
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Table 3. Comparative analysis of MHP potential site 


Description Unit CENENE- i o ee ee ee 
Estimated Measured Designed Estimated Measured Designed 
Catchement Area ha 1,595 6,515 
Potential Flow m?/s 1.30 1.21 0.80 10.20 8.88 6.00 
Potential Head m 25.00 22.00 20.60 125.00 120.00 110.00 
Potential Power kW 255.00 182.00 129.00 10,000.00 8,365.00 5,180.00 


The differences between the hydropower estimation and measurement are 39% and 20% for Site 1 
and Site 2 consecutively. The proposed hydropower design considers the minimum flow and risk factors for 
changes that may occur. There fore the head and flow rate design are set to be smaller than that of the 
estimation and measurement. 

Remote sensing interpretation generally produces overestimated values compared to the 
measurement data due to less detail of the data resolution and the map scale used. However, the results of the 
remote sensing interpretation can give an idea of the existing potential so that it so can be used as a basis for 
assessing the hydropower potential in an area, especially the ones with inadequate measurement data. 
Detailed measurements in the field are needed to validate data and obtain more precise results. 


4. CONCLUSION 

This paper discusses the potential study of MHP with the help of Landsat 8 OLI satellite data. 
Sukaati watershed, Cibuni river, West Java was chosen because it has great potential for electrical energy 
used for industry or society. The methodology integrates interpretation of satellite imagery and field surveys 
related to geological conditions, topography, land cover, environment, and geological disasters such as 
landslides. The main constituent rocks are volcanic rocks, namely breccia and andesite lava, with land cover 
dominated by tea plantations. The low potential for landslides greatly supports the construction of the MHP. 
The results of satellite data interpretation of geological conditions, land cover, topography, and geological 
disasters show two potential locations for MHP, suitable for hydroelectric power plants with 129 kW and 
5.18 MW hydropower designs. However, detailed surveys and field measurements are still needed to validate 
the interpretation and analysis of remote sensing imagery results. 
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